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ABSTRACT 

We investigate the cold and warm dust properties during galaxy interactions using 
a merging galaxy sample ordered into a chronological sequence from pre- to post- 
mergers. Our sample comprises a total of 29 merging systems selected to have far- 
infrared and sub-millimeter observations. The sub-millimeter data are mainly culled 
from the literature while for 5 galaxies (NGC3597, NGC3690, NGC6090, NGC6670 
and NGC 7252) the sub-millimeter observations are presented here for the first time. 
We use the 100~to-850/.im flux density ratio, /looZ/sso, as a proxy to the mass fraction 
of the warm and the cold dust in these systems. We find evidence for an increase in 
/ioo//850 along the merging sequence from early to advanced mergers and interpret 
this trend as an increase of the warm relative to the cold dust mass. We argue that 
the two key parameters affecting the /looZ/sso flux ratio is the star-formation rate 
and the dust content of individual systems relative to the stars. Using a sophisticated 
model for the absorption and re-emission of the stellar UV radiation by dust we show 
that these parameters can indeed explain both the increase and the observed scatter 
in the /looZ/sso along the merging galaxy sequence. We also discuss our results under 
the hypothesis that elliptical galaxies are formed via disc galaxy mergers. 
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1 INTRODUCTION 

Early studies on the dust properties of spiral galaxies us- 
ing IRAS data suggested a gas-to-dust ratio about 5-10 
times larger than the Milky Way (Devereux & Young 1990; 
Sanders et al. 1991). This apparent conflict was interpreted 
as indirect evidence for the presence of a cold dust com- 
ponent at a temperature T=15-25K heated by the diffuse 
interstellar radiation field (Cox, Krugel & Mezger 1986; 
Rowan- Robinson & Crawford 1989). Dust at such low tem- 
peratures remains undetected at the IRAS wavebands, that 
primarily probe warmer dust (T > 30 K), and can only be 
identified by observations at longer wavelengths (> lOO/im). 
It was only recently that the SCUBA bolometer array on 
the James Clerk Maxwell Telescope (JCMT) has provided 
sufficient sensitivity at submillimeter (sub-mm) wavelengths 
(450, 850/im) to directly confirm the presence of a domi- 
nant cold dust component in nearby spirals (Alton et al. 
1998a; Alton et al. 1998b; Frayer et al. 1999; Papadopoulos 
& Seaquist 1999; Dunne et al. 2000; Dunne & Bales 2001). 

In more active luminous infrared galaxies (LICs), with 
enhanced star- formation or the presence of a central AGN, 
sub-mm observations also suggest the presence of cold cirrus. 
Unlike quiescent spirals however, LIGs also show a promi- 



nent warm dust component (T=30-60 K) most likely heated 
by the starburst (Klaas et al. 2001; Lisenfeld, Issak & Hills 
2000). The enhanced activity in LIGs is believed to be trig- 
gered by interactions and mergers. It is therefore, interesting 
to explore the role of such violent events in heating the dust 
and hence, modifying the relative amounts of the cold and 
the warm dust component in galaxies. Klaas et al. (2001) 
have searched for systematic differences in the far-infrared 
(far-IR) to sub-mm luminosity ratio of interacting and non- 
interacting LIGs classified on the basis of their optical mor- 
phology. Their analysis did not show any trends suggest- 
ing either problems with their classification scheme or sim- 
ilar dust temperature distributions in the two sub-samples. 
Studies of selected interacting galaxies (mostly LIGs) in- 
deed, suggest that cold dust is present in these systems even 
in the regions experiencing starburst activity where the stel- 
lar radiation field heating the dust is most intense (Haas et 
al. 2000). This may be due to effective shielding of the cold 
dust component by warmer dust. 

In this paper we further explore the issue of dust heat- 
ing during galaxy encounters by exploiting the large SCUBA 
database as well as sub-mm observations presented in the 
literature. Our interacting galaxy sample is compiled from 
the Georgakakis et al. (2000) study and comprises well sepa- 
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rated spirals, systems close to nuclear coalescence and young 
merger remnants. Our sample is larger compared to previous 
studies and also spans a wider range of interactions stages. 

Section 2 discusses our sample selection while Section 3 
presents the reduction of the data used in our study. In Sec- 
tion 4 we present our results which are discussed in Section 
5. Finally section 6 summarises our conclusions. Throughout 
this paper we adopt Ho = 75kms~^ Mpc~^. 



2 THE SAMPLE 

In this paper we use the interacting galaxy sample presented 
by Georgakakis et al. (2000) comprising disk galaxy mergers 
of similar mass. This is culled from a number of studies (Keel 
& Wu 1995; Gao & Solomo 1999; Gao et al. 1998; Surace et 
al. 1993) in an attempt to minimize any biases introduced by 
different selection criteria (e.g. far-IR, morphological selec- 
tion) used to compile interacting galaxy samples. We note 
however, that most of the galaxies in that study are far- 
IR luminous and are also biased against mergers occurring 
along our line of sight. Therefore, although the Georgakakis 
et al. (2000) sample is by no moans statistically complete, 
it could be regarded as representative of interacting systems 
and merger remnants spanning a wide range of properties. 
The sample employed in this study is presented in Table 
1 and comprises only those systems from the Georgakakis 
et al. (2000) study that have sub-mm data from either the 
SCUBA archive or the literature. 

Following Georgakakis et al. (2000) we use the galaxy 
'age' parameter providing an estimate of the evolutionary 
stage of each system relative to the time of the merging 
of the two nuclei. Such a chronological sequence allow us 
to explore the evolution of the galaxy properties at dif- 
ferent times during the interaction. Negative 'ages' are for 
pre-mergers while positive 'ages' correspond to merger rem- 
nants. For pre-mergers the 'age' is estimated by dividing 
the projected separation of the two imclci, Sr, by an (arbi- 
trary) orbital decay velocity v = 30kms~'^. It is clear that 
the 'age' parameter for pre-mergers is affected by projection 
effects or different interaction geometries. However, to the 
first approximation, it provides an estimate of the stage of 
the merging and allows plotting of pre- and post-mergers 
on the same scale. For post-mergers we adopt the evolution- 
ary sequence defined by Keel & Wu (1995) using dynamical 
and morphological criteria. In particular, the 'age' param- 
eter for these systems is calculated by multiplying the dy- 
namical stage number, defined by Keel & Wu (1995), by the 
factor 4 x 10® yr. This conversion factor is found to be appro- 
priate for the 3 merger remnants in the Keel & Wu sample 
with available spectroscopic estimates (i.e. NGC 2865, NGC 
3921, NGC 7252; Forbes, Ponman & Brown 1998). It should 
be stressed that the 'age' parameter for both pre- and post- 
mergers does not represent an absolute galaxy age but is an 
indicator of the evolutionary stage of the interaction. 



3 OBSERVATIONS 

In addition to the sub-mm data available in the litera- 
ture we have analysed unpublished SCUBA archival data 
for five merging systems in the present sample: NGC 3597, 



NGC 3690, NGC 6090, NGC 6670 and NGC 7252. The obser- 
vations of these galaxies were made with the SCUBA sub- 
mm bolometer array (Holland et al. 1999) which is mounted 
on the JCMT and provides simultaneous imaging at 450 
and 850 /im for a region of sky of 2.3arcmin in diameter. 
SCUBA is an array of 37 bolometers (HPBW=14.7arcsec) 
at 850 fj,m and 91 (HPBW=7.5 arcmin) at 450 //m. In order 
to provide fully sampled images, the secondaxy mirror moves 
in a 64-step jiggle pattern, with the integration time lasting 
1 s at each position. After the 16 stops of jiggle pattern, the 
telescope nods in order to allow for slowly varying sky gra- 
dients. For NGC 3597 (observed in 5- Apr- 2000), NGC 3690 
(observed in 24-Nov-98), NGC 6670 (observed in 5, 6-Nov- 
97) and NGC 7252 (observed in 5-Jul-97) the chop throw 
was 120, 80, 80 and 120 arcsec respectively. For NGC 6090 a 
chop throw of 120, 80 and 60 arcsec was used during the ob- 
serving runs of 12-Mar-1998, 5-Dec-1997 and 4-Jan-2000 re- 
spectively. All observations were made with grade- 1 weather 
(r225 < 0.05) except for 5-Apr-2000 and 4-Jan-2000 where 
weather conditions were of grade-3 {T225 ~ 0.1). 

The data were reduced using the SURF software pack- 
age (Jenness et al. 1997). We made corrections for the atmo- 
spheric absorption using the opacities derived from skydip 
measurements taken regularly every night. Noisy bolometers 
were flagged and large spikes were removed. During each 
night, several calibration sources were mapped in order to 
flux calibrate our images. The SCUBA beams have moder- 
ate error lobes and since all our targets arc extended sources, 
we corrected the calibration in order to account of emission 
outside the central beam. We estimated the flux loss due 
to the secondary lobes by evaluating the average peak-to- 
aperture flux ratios measured for the calibration sources for 
each observing night. All five new galaxy systems presented 
in this work were detected at 850 jim while only NGC 3690 
showed significant signal in the 450 /um band. For all our ob- 
servations wc estimated a calibration uncertainty (based on 
differences among the nightly calibration values) which are 
reported in Table 1. The 850 pm emission of these galaxies 
is presented in Figure 1 with contours overlaid on top of 
the optical DSS images. Unfortunately, the relatively large 
beam size of SCUBA (~ 4 kpc at a distance of 50 Mpc) does 
not allow for a morphological study of these systems. We 
note however, that some of these systems (mostly early pre- 
mergers) have 850 /xm emission that shows two peaks, while 
others (the more evolved systems) show a single peak. 



4 RESULTS 

The dust properties of our sample are quantified using the 
100-to 850 //m flux density ratio, /looZ/sso- This provides 
an estimate of the relative mass fraction of the warm and the 
cold dust component. Indeed, the sub-mm wavelengths (e.g. 
450 and 850 pm) are sensitive to cold dust (T « 15 — 25 K; 
Alton et al. 1998a; Alton et al. 1998b; Bianchi et al. 1998; 
Dunne et al. 2000; Dunne & Bales 2001) while the 100 /um 
flux density primarily probes the warmer dust (T « 35 K) 
heated by the star-formation activity (Devereux & Young 
1990; Devereux et al. 1994; Klaas et al. 2001; Misiriotis 
et al. 2004). Also the presence of an AGN in some of our 
systems does not affect our results. This is because at the 
wavelengths used in this study the contamination of the ob- 
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Figure 1. Sub-mm images of NGC 3597, NGC 3690, NGC 6090, NGC 6670 and NGC 7252. Contour levels are 0.020, 0.033, 0.047, 
0.060, 0.073 Jy/bearn for NGC 3597; 0.030, 0.068, 0.107, 0.146, 0.184, 0.222 Jy/bearn for NGC 3690; 0.020, 0.051, 0.082 Jy/beam for 
NGC 6090; 0.040, 0.063, 0.086 Jy/beam for NGC 6670 and 0.010, 0.017, 0.025, 0.032, 0.040 Jy/beam for NGC 7252. 



served fluxes due to AGN emission is expected to be small 
(e.g. Delia Ccea et al. 2002; Zezas, Ward & Murray 2003). 
This is also discussed below. 

We further explore the connection between /ioo//s50 
and the dust content by estimating the warm (Mw) and 
the cold (Mc) dust mass of our sample galaxies and then 
plotting /looZ/sso against the ratio Mc/Mw- The Mw, Mc 
are calculated by fitting two modified Planck functions, each 
with different temperature, to the observed far-IR to sub- 
mm spectral energy distribution of our sources. We follow 
the method described by Dunne & Bales (2001) keeping the 
dust emissivity index fixed to /3 = 2. In this exercise we 
combine the IRAS and the SCUBA data presented in Table 1 
with the flux density measurements at other far-IR and sub- 
mm wavelengths from the literature. Only systems with data 
in at least four wavelengths, of which at least two should lie 
in the sub-mm regime, are fitted. 

The emission at a particular wavelength is modeled by 
a two-temperature graybody described by 

Fx = ^B{\,Tw) + ^B{\,Tc), (!) 

where Tw and Tc are the temperatures of the warm and 
the cold dust component respectively and /3 is the dust 
emissivity index, fixed to /3 = 2 for consistency with the 
study of Dunne & Bales 2001. Nw and Nc are normaliza- 
tion constants that determine the mass of each of the warm 
and cold components respectively. As already argued in Al- 
ton et al. (1998a), although the dust masses that are cal- 
culated with this technique are not very accurate (due to 



the adopted dust properties), the cold-to-warm dust mass 
ratio Mc/Mw (which is equivalent to Nc/Nw) provides a 
more secure estimate of the relative content of warm and 
cold dust within the galaxies. We also estimate the statisti- 
cal error of the Mc /Mw ratio due to the uncertainty in the 
measured fluxes using the bootstrap method described by 
Dune & Bales (2001). For each galaxy in our sample mock 
SBDs were constructed by randomising the fluxes at each 
wavelength, assuming that they follow a Gaussian distribu- 
tion with mean the measured flux and standard deviation 
the uncertainty estimated from the observations. A total of 
100 mock SBDs were produced for individual galaxies in 
our sample and each of them was then fltted with two mod- 
ified black bodies in the same way as with the real data. 
The standard deviation of Mc/Mw ratio from the 100 trials 
provides an estimate of the uncertainty in this parameter. 
We note however, that this uncertainty does include system- 
atic biases that may affect the Mc /Mw ratio. For example, 
the sparse wavelength coverage, uncertainties in the opacity 
coefficients Kd and the emissivity index /3 as well as degen- 
eracies in the fitting parameters make the decomposition of 
the observed SBD into two blackbodies not unique (Klaas 
et al. 2001). 

The results of the fits are presented in Table 2. Figure 
2 plots /ioo//85() against the ratio Mc/Mw for the merg- 
ing galaxies presented in this study (solid circles) and the 
galaxies presented in Dunne & Bales 2001 (open circles). De- 
spite the large scatter there is evidence for a correlation in 
the sense that sources with a lower cold-to-warm dust mass 
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Table 1. Observational information of the galaxies. 
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a: Chini et al. 1986; b: Krucgcl ct al. 1988; c: Ealcs, Wynn- Williams & Duncan 1989; d: Carico ct al. 1992; 

e: Clements, Andrcani & Chase 1993; f: Rigopoulou, Lawrence & Rowan- Robinson 1996; g: Benford et al. 1999; h; Prayer et al. 1999; 

i: Haas et al. 2000; j: Dunne & Ealcs 2001; k: Klaas ct al. 2001; 1: this work.; m: Calzetti et al. 2000; n: Spinoglio, Andreani & Malkan 2002; 

o; Andreani & Pranceschini 1996; p; Siebenmorgen, Krugel & Chini 1999 



fum/.fsr,() and the age parameter arc uncorrelated can be 
rejected at the « 3.5a confidence level. 

We note that this statistically significant correlation is 
not due to /loo or fsso being individually correlated with 
the age parameter. We normalise /loo and /sso with the K- 
band flux {/k) of our systems from the Two Micron All Sky 
Survey (2MASS; .Jarrett et al. 2000). Using the Spearman 
rank correlation test we find that the probability that the 
age parameter is uncorrelated with /ioo//k and /s5o//if of 
4 (w 2a) and 28 (« la) per cent respectively. Therefore, the 
null hypothesis that there is no correlation between these 
parameters cannot be rejected at a high confidence level. 
Also, the trend between /looZ/sso and the age parameter 
is not affected by the presence of an AGN in some of our 
sample galaxies. This is demonstrated in Figure 3 where 
systems identified as AGN in the catalogue of Veron-Getty & 
Vcron (2003) arc plotted with crossed circles. The AGNs in 
that figure follow a similar trend to the remaining systems. 

The mean /100/./850 ratio of isolated spirals from the 
sample of Misiriotis ct al. (2004) is also plotted in Figure 3 
at an arbitrary age of -13. Based on the discussion above 
this trend can be interpreted as an increase of the warm rel- 
ative to the cold dust mass from isolated spirals and early 



fraction (low Mc /Mw values) also have enhanced /loo //sso 
fiux density ratio. The Spearman's rank correlation test re- 
jects the null hypothesis that the /looZ/sso and the Mc /Mw 
are uncorrelated with a probability 0.015 corresponding to 
« 2.5a. Similar conclusions are obtained if we use the larger 
sample of Dunne & Eales (2001) that comprise both inter- 
acting systems and isolated galaxies. 

The large scatter in Figure 2 may be intrinsic suggest- 
ing that parameters other than Mc/Mw may affect the 
/ioo//850. However, the systematics discussed above may 
also bias the Mc /Mw ratio and contribute to the observed 
scatter. We avoid the above systematics by using observ- 
ables to describe the dust properties (e.g. /loo/Zaso) rather 
than parameters derived from spectral fittings. 

Figure 3 plots the ratio /looZ/sso as a function of the 
age parameter. As explained in Section 2, negative ages cor- 
respond to pre-merger stages, zero is for the nuclear coales- 
cence while positive ages are for merger remnants. In Fig- 
ure 3 there is evidence for an increase in the /loo //sso ra- 
tio as the interaction progresses toward the final stages of 
the merger. We estimate a Spearman rank correlation coeffi- 
cient r = 0.64 and a probability that there is no correlation 
P = 0.07 per cent. Therefore, the null hypothesis that the 
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Table 2. Dust temperatures and cold-to-warm dust mass ratio 
derived by fitting a two-temperature graybody curve with /3 = 2. 
The uncertainty of the cold-to-waxm dust mass ratio is given in 
tha last coloumn. 
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interacting systems to advanced mergers. Such a transition 
of the cold dust to warmer temperatures can be attributed 
to heating by the more intense UV radiation field of the 
enhanced star-formation triggered by the merging (e.g. Ca- 
soh et al. 1991; Keel & Wu 1995; Georgakakis et al. 2000). 
However, the lack of a statistically significant correlation 
between /ioo//if (« 2cr; providing an estimate of the SFR 
normalised to the total galaxy mass) and age suggests that 
parameters other than the SFR also play an important role 
in moderating the /looZ/sso ratio. 

We further explore this using the sophisticated model of 
Misiriotis et al. (2004) for the absorption of the UV/optical 
radiation by dust and its re-emission at far-IR and sub- 
mm wavelengths. This model produces UV to sub-mm SEDs 
for a given input system geometry, stellar population (e.g. 
young and old stars), dust content and total galaxy mass. 
Our goal is not to derive physical parameters for the individ- 
ual interacting systems but to show that this model, under 
simple assumptions, can reproduce the observed /looZ/sBO 
range of our sample. 

The interacting galaxies studied here clearly have com- 
plex morphology that is difficult to quantify and to model. 
To the first approximation we assume simple exponential 
disk-galaxy geometry. Although clearly a more complex ge- 
ometry is required for each one of our interacting systems, 
this is an approximation that allows some insight into how 
different parameters (e.g. star-formation rate, dust content) 
affect the /loo / /sso ratio. Since wo do not attempt to model 
individual systems the exact geometry of the interaction is 
not a crucial parameter for this study. The total mass of the 
model galaxy is parametrised by its i^T-band luminosity. We 
adopt Mk ~ —24.5 mag, the mean absolute magnitude of 
our sample galaxies. The dust mass of the fiducial galaxy is 
assumed to be in the range w lO'^ — 10^ M0 , similar to the 
dust masses estimated using the simple black-body fits de- 



scribed above. The stellar content is directly associated with 
the density of the UV radiation responsible for the heating 
of the dust and is parametrised by the star-formation rate. 
Here we use three different values for the star-formation rate 
SFR=1, 10 and 60 Mq- For a given SFR and dust mass (in 
the range « 10^ — 10** Mq) we derive the /looZ/sso ratio pre- 
dicted by the Misiriotis et al. (2004) model. The shaded re- 
gions in Figure 3 represent the /loo //sso for different SFRs. 
The lower and upper boundaries of these regions correspond 
to dust masses of 10^ and 10* Mq respectively. Clearly, both 
the SFR and the total dust mass (relative to the stellar mass) 
are important for the /loo //sso ratio. 

As expected the model predicts that higher star- 
formation rates heat the dust resulting in higher /ioo/Zbso 
ratios. Indeed, the model of Misiriotis et al. (2004) predicts 
an increase in the mean (weighted by mass) dust temper- 
ature of the cold component for different SFRs: we find 
that for SFR=1 Mq/jt the dust temperature is between 16 
and 18 K, for WMq/jj: it is between 19 and 22 K while for 
SFR=6OM0/yr it is between 24 and 28 K. 

Also, keeping the SFR fixed, more dusty systems have 
lower /100//8SO which can be attributed to higher fraction 
of UV/optical photons absorbed by the dust, heating it to 
temperatures T ^ 30K and then being re-emitted primar- 
ily at sub-mm wavelengths (e.g. SSO/nm). For example in 
the case of the Misiriotis et al. (2004) study, the fraction 
of UV/optical photons absorbed by the diffuse dust com- 
ponent in individual systems in their sample (as predicted 
by their best fit model) is anti-correlated with the observed 
/loo/Zaso at the 2.5cr confidence level. This may also partly 
explain the scatter seen in Figure 3 since the relative fraction 
of the warm-to-cold dust for a given SFR varies substan- 
tially with dust content. Other factors such as projection 
effects, differences in the details of individual interactions 
(e.g. geometry, initial conditions, bulge-to-disk ratio) may 
also contribute to the scatter in Figure 3. We also note that 
the diffuse UV/optical radiation field does not have a large 
impact on the lOO^m emission, which is primarily domi- 
nated by starburst activity taking place in localised galaxy 
regions. 

At post-merger stages the evolution of /100//850 with 
age after the merging event remains poorly constrained since 
only 3 systems have SCUBA measurements available. This 
is also a problem for evolved ellipticals, the alleged product 
of disc galaxy mergers, with dust properties (e.g. tempera- 
ture, content) that remain largely unknown. This is mainly 
because these systems are believed to have very little if any 
dust and hence are difficult to detect at far-IR and sub-mm 
wavelengths. Previous observations of ellipticals at sub-mm 
wavelengths targeted far-IR luminous and optically pecu- 
liar systems to maximize the probability of detection (e.g. 
Fich & Hodge 1993; Wikhnd & Henkel 1995; Leeuw et al. 
2004). As a result the /ioo//sso ratios of these systems may 
not be representative of ellipticals. We attempt to constrain 
the /ioo//s5o ratio of evolved ellipticals using the sample of 
Temi et al. (2003). These authors use ISO observations to 
constrain the far-IR (60-200 /im) SEDs of early type galax- 
ies (not selected on the basis of their far-IR luminosity) and 
to estimate their warm and cold dust content using meth- 
ods similar to those outlined here (e.g. modified black-body 
fittings) . Although their results are somewhat uncertain due 
to the lack of sub-mm observations which are critical as a 
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Figure 2. fum/fsso flux ratio against the cold-to-warm dust 
content. Tlie filled circles are the interacting systems in the 
present sample while the open circles are the galaxies presented 
in Dunne & Bales (2001). 



probe of the cold dust, this study provides the only con- 
straints on the dust properties of early type systems using 
data extending beyond lOO^im. Using the best fit model of 
Temi et al. (2003) we extrapolate their observed SEDs to 
850 fjxn and estimate the mean /loo //sso flux density ratio 
for the elliptical galaxies in their sample. This mean value 
(54 ±33) is plotted in Figure 3 (the star arbitrary plotted at 
an age of +5) and should be considered with caution since 
the /sso flux density is not directly observed. Nevertheless, 
taken at face value it suggests that the ratio /ioo/./85() in 
evolved ellipticals is about 1 dex lower compared to young 
merger remnants. This can be interpreted as dust cooling 
after the merger event under the hypothesis that elliptical 
galaxies are formed via disc galaxy mergers. We caution the 
reader however, that other mechanisms such as spattering 
of dust grains by X-ray radiation as well as the dispersal of 
dust in star formation may also affect the /looZ/sso ratio. 
Observations at sub-mm wavelengths of elliptical galaxies 
and post-mergers are essential to further explore different 
scenarios. 



5 DISCUSSION 

In this paper we use an interacting galaxy sample ordered 
into a chronological merging sequence to explore dust prop- 
erties during gravitational galaxy encounters. We find evi- 
dence for an increase in the flux density ratio /loo //sso from 
isolated spirals and early interacting systems to intermedi- 
ate and late stage mergers close to nuclear coalescence. We 
interpreted this trend as an increase of the warm relative 
to the cold dust mass. Both the star-formation triggered by 
the merging and the dust content of individual systems are 
important parameters for modifying the fum / /sso ratio and 
the observed trend along the merging galaxy sequence. 

Klaas et al. (2001) have also investigated the far-IR to 
sub-mm properties of luminous infrared galaxies and have 
searched for systematic differences between isolated and in- 
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Figure 3. /loo/Zsso flux ratio against the age parameter for our 
sample of interacting galaxies. The filled circles are the interact- 
ing galaxy sample in this study. A cross on top of symbol sig- 
nifies an AGN in the catalogue of Vcron-Cctty & Vcron (2003). 
The square corresponds to the mean /looZ/sso of isolated spirals 
from the sample of Misiriotis et al. (2004) arbitrarily plotted at 
an age parameter of -13. The star is the mean /ioo//sso of ellip- 
ticals from the sample of Temi et al. (2003; see text for details) 
arbitrarily plotted at an age parameter of -|-5. 



teracting systems classified on the basis of their optical mor- 
phology. Unlike our results they find no differences in the 
relative amounts of warm and cold dust between the two 
samples. However, these authors do not attempt to order 
their interacting systems on the basis of the angular sepa- 
ration between the two nuclei (ie. the age parameter). Also, 
their sample spans a small range in bolometric luminosity 
and therefore well separated early interacting systems may 
be underrepresented. 

Dunne & Eales (2001) have performed the most exten- 
sive study of the sub-mm dust properties of nearby galaxies 
todate using the SCUBA bolometer array at the JCMT. 
They find that in quiescent spirals most of the warm dust 
(T > 30 K) is concentrated in nuclear regions with the cold 
dust component (T = 15 — 20 K) dominating the outer disc. 
The situation is strikingly different in starbursts where they 
find that the warm dust region is more extended and there- 
fore dominates the SED. This is attributed to the more 
intense interstellar radiation field in actively star-forming 
galaxies. Although Dunne & Eales (2001) do not discuss 
possible differences in the sub-mm properties between inter- 
acting and isolated galaxies in their sample, most of our in- 
teracting/merging systems, particularly the more advanced 
ones, are in a starburst phase. Therefore, on the basis of 
the Dunne & Eales (2001) findings, for these galaxies we 
would expect a significant warm dust component compared 
to isolated (on average more quiescent) spirals. 

We note however, that the evidence above does not 
suggest the absence of cold dust in interacting starbursts. 
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Indeed, detailed study of selected systems (e.g. ARP 220, 
NGC 4038/4039) shows significant amounts of cold dust that 
may dominate the total dust content (Klaas et al. 1997; Hcias 
et al. 2000: Klaas et al. 2001; Dumkc, Krause & Wielebinski 
2003). Dunne & Ealcs (2001) argue that the most active star- 
burst activity is usually taking place in well-shielded dust 
enshrouded regions allowing the presence of very cold dust 
at the same regions. 

Moreover, using our interacting galaxy sample we at- 
tempt to explore possible connections between the dust 
properties of isolated spirals, mergers and evolved ellipticals. 
This is particularly important in the light of the merger hy- 
pothesis postulating that ellipticals are the product of disc 
galaxy mergers. However, our sample comprises only a few 
post merger remnants while the cold and warm dust prop- 
erties of evolved ellipticals remain poorly constrained. Nev- 
ertheless, the tentative result of our study is that young 
merger remnants have /ioo//s50 ratios about 1 dex higher 
than evolved ellipticals. This may suggest that the dust cools 
with time after nuclear coalescence, although other scenar- 
ios can also explain the lower /ioo//850 ratios in these sys- 
tems. Future observations using sensitive instruments such 
as ALMA and the Spitzer satellite will allow us to probe the 
largely unexplored far-IR to sub-mm regime of early type 
galaxies. These data have the potential to accurately char- 
acterize, for the first time, the dust properties of ellipticals 
and merger remnants and to search for possible evolutionary 
eflfects between these systems. 



6 CONCLUSIONS 

In this paper we explore the cold and warm dust properties 
of gala^xies during gravitational encounters using a merging 
galaxy sample ordered into a time sequence and compris- 
ing well separated spirals, systems close to nuclear coales- 
cence and young merger remnants. We use the 100 to 850/im 
flux density ratio, /ioo//s50, as a proxy to the mass fraction 
of the warm and the cold dust in these systems. We find 
evidence for an increase in /loo/Zsso along the merging se- 
quence from isolated spirals and early interacting systems to 
advanced mergers. We interpret this trend as an increase of 
the warm relative to the cold dust mass and argue that the 
star-formation and the dust content of individual systems 
are important for modifying the /loo //sso ratio. 
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